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■"■•^ ' Our quasi-particle model for deconfined matter near T c is reviewed. The extrap- 

' olation of lattice QCD data to a finite baryo-chemical potential is discussed. De- 

Oh, termined by the chiral transition temperature T c , the resulting equation of state 

of neutral and /3-stable deconfined matter is soft and limits size and mass of pure 
quark stars. 
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1. Introduction 



Numerical Monte Carlo evaluations of QCD have succeeded to deliver the 
equation of state (EoS) of strongly interacting matter (quarks and gluons) 
with increasing accuracy. Near and slightly above the chiral transition 
temperature (which seems to coincide with the deconfinement temperature) 
T c , the EoS shows a non-trivial behavior. It is, therefore, a challenge to 
develop models based on QCD, and catching its relevant degrees of freedom, 
which describe the EoS in a quantitative and transparent manner. A big 
deal of progress has been made by starting from first principles and arriving, 
after a chain of systematic approximations, at models which describe the 
lattice QCD data. Despite of several attempts, the range of applicability 
is restricted to temperatures T > 2.5T C . 1,2 ' 3 Here, we are interested in 
the range T c — 3T C which covers the physically important range for heavy- 
ion physics. We are going to review our phenomenological quasi-particle 
model and to describe its use for mapping lattice QCD data to a finite 
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baryo-chcmical potential. In contrast to the models derived directly from 
QCD, 1:2,3 we use lattice QCD data to adjust a few parameters. As an 
application we discuss implications for cold pure quark stars. 



2. Quasi-particle model 

Our model 4 ' 5 is based on the entropy density as combinatoric quantity 
measuring the population of states, s = J2i= q g s id ( m » M) > T, /i) , which 
implies for the pressure p — J2i= q g Pid ~ B(T, /i) via s = dp/dT. The sub- 
scripts "id" indicate the corresponding ideal gas expressions. All interac- 
tions are assumed to be included in the quasi-particle masses m 2 (T, /i) oc a s , 
which correspond to asymptotic masses of excitations near the light cone 
with explicitly given T and /i dependence from HTL/HDL calculations, 
and the effective coupling strength a s (T, /i) is thought to include non- 
perturbative effects beyond that. Expanding p(T), e.g., for the pure gluon 
plasma with the lowest order rrii = ^<7 2 T 2 one recovers the known perturba- 



tive series, p = psb (l 
= 127r/(33-2AT / )lnf 
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•) . With the parameterization 



we can describe available lattice QCD at 
/i = 0. Fig. 1 exhibits two examples. 

Due to thermodynamic consistency the trace anomaly is also described 
correctly. Opposed to the quasi-particle masses, the screening mass drops 
near T c . 4 The analysis 6 of the temporal quark and gluon propagators sup- 
ports our model. 
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Figure 1. Comparison of lattice QCD results and our quasi-particle model. Left: pure 
gluon gas, data from Ref. 7 ; Right: Nf = 2 quark-gluon plasma, data from Ref. 8 ; open 
(closed) symbols are for light (heavier) quarks; the shadowed region depicts a continuum 
estimate 9 ; the hatched area indicates the pure gluon plasma from the left panel. The 
parameters T s , A are adjusted to lattice QCD data. For further details see Ref. J . 
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Other quasi-particle formulations are conceivable (see Refs. 4 ' 10 for quo- 
tations). For example, in Ref. 11 a different handling of the quasi-particle 
masses has been advocated, where the degeneracies are taken as tempera- 
ture dependent and adjustable quantities. Our formulation is based on the 
(/>-derivable formalism, where the entropy density is a particularly interest- 
ing quantity to be considered. 12 

Ref. 13 makes a strong claim that the structure of the above quasi- 
particle formulation is not conform to the exactly solvable QCD in the 
large Nf limit; however, as argued in Ref. 14 this disclaimer may not be 
relevant for the present model which is not aimed at covering this limit of 
QCD. 



3. Extrapolation to finite baryo-chemical potential 

The above discussion is constrained to vanishing baryo-chemical poten- 
tial. Due to the stationarity condition of the thermodynamical potential, 
dp/drrij = 0, the entropy and baryon density are given by the sum of the 
individual quasi-particle contributions, s, = dpi/dT and rij = dpi/d^. On 
the other hand, the pressure p must be a thermodynamic potential so that 
ds/d/i — dn/dT. This results in a partial differential equation for d(T, /n), 
which can be solved as well-posed boundary problem with a(T,/i — 0) ad- 
justed to lattice QCD data. A solution by the method of characteristic 
curves is exhibited in Fig. 2 together with the resulting pressure. 




12 3 4 



fj,/ T c 

Figure 2. Left: The characteristic curves in the T — fj, plane. The heavy dashed curve 
is the estimated 15 phase border line (uncertainty indicated by the hatched area). Left 
to the dash-dotted curve the pressure becomes negative. Right: The EoS in our quasi- 
particle model (T and fi in units of T c , p in units of psb)- Further details in Ref. 5 . 
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Very recently, lattice QCD calculations with // > became available. 16 
This offers the chance to check the validity of our extrapolation of the EoS. 
A first result is the curvature of the critical line at T c , 15,16 see Fig. 2. Some 
idea on the robustness of our extrapolation can be gained by a comparison 
to the HTL quasi-particle model applying the same method. The results 
in Ref. 17 show indeed only a small numerical deviation from our results. 



4. Pure quark stars 

Soon after the discovery of quarks and gluons as the constituents of hadrons, 
the existence of quark cores in neutron stars or pure quark stars has been 
conjectured. Since that time a lot of work has been done to explore various 
consequences of this hypothesis. If the star matter exhibits a first-order 
phase transition above nuclear density, a third family of cold quark stars 
can appear beyond the stable branches of white dwarfs and ordinary neu- 
tron stars. 18 The transition, triggered e.g. by accretion, of a neutron star to 
a quark star can be accompanied by an ejection of the outermost layers. 19 
This may result in particular supernova events. 20,21 Recently, the possible 
occurrence of the third family, which would allow so-called twin star config- 
urations, has been rediscovered. 22 The necessary discontinuity in the EoS 
may be realized by the chiral transition, but also other phase changes of 
nuclear matter may allow twin stars. 21 

From the universality of the EoS at fi = 0, the extrapolated EoS at 
fi > is rather robust, supposed there is no change of the quasi-particle 
structure. All our examined examples deliver approximately the relation 
e = 4B + ap with a = 3 • ■ ■ 4 and B 1 / 4 = (3T C , j3 = 1.7 ■ ■ ■ 2.1 at T = 0. 
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Figure 3. Mass-radius relation of cold 
pure quark stars delivered by our EoS. 



Integration of the TOV equations 
results in light and small quark 
stars, 5 see Fig. 3. Such objects are 
not (yet?) found. But one should 
have in mind, that the outermost 
layers of our quark star models are 
unstable with respect to a phase 
conversion. Continuing the EoS in 
the region where the pressure of our 
model vanishes by another EoS with 
larger pressure can change signifi- 
cantly the mass-radius relation. 
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5. Summary 

In summary we have reviewed our quasi-particle model of deconfmed mat- 
ter. To test in more detail the reliability of our model we still need more 
information on the EoS for several flavor numbers with systematic con- 
tinuum extrapolation and with realistic quark masses. Nevertheless, our 
model allows already at the present stage an extrapolation of the EoS to 
cold quark matter. The resulting cold pure quark stars are light and small. 
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